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ABSTRACT 


This thesis discusses the conversion of the MOD6DF computer program 
for use on the IEM-360 computer at the Naval Postgraduate School. The 
functioning program was modified to investigate the impact miss distance 
for the Supersonic Tactical Missile. When the initial y-displacement 
error exceeded 1800 feet, the missile did not acquire the target. All 
errors smaller than this resulted in miss distances within 0.5 feet of 
the target. The midcourse guidance reference altitude was changed to 
reflect a sea-skimming missile. This simulation ran and impact was 
recorded. An attempt at adding random noise to the homing seeker was 
tried, but revealed that more information is required on this topic. 


The MOD6DF computer program was successfully converted and altered to 


run using the simplified ramjet model. 
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TABLE OF SYMBOLS AND ABERSVIATIONS 


The following is a list of the abbreviations and some of the more 
common fortran symbols used in the MOD6DF computer program. Each symbol 
is defined in two ways. The primary source of identification is by its 
COMMON( 3415) location and the second is by its fortran symbol. The 
fortran symbol is not always a good identifier since it will change from 
subroutine to subroutine (i.e. W and WE both are used for missile 


weight ) e 


A. ABEREVIATIONS 


Symbol Definition 

LP tangent plane axes 
RBA body axes 

SA stability axes 


B. FORTRAN SYMBOLS 


Symbol Definition 

TXBA(073) missile thrust in the x-direction in BA (1b) 
TYBA (074) missile thrust in the y-direction in BA (1b) 
TZBA(075) missile thrust in the z-direction in BA (1b) 
W( 086) missile weight (1b) 

$(110) missile reference area (ft*) 

CA(111) drag coefficient 

CY( 113) side-force coefficient 

CZ(115) normal force coeffieient 

CBAR(116) mean aerodynamic chord (ft) 

CMQ(118) damping in pitch coeffieient 

CNR(119) damping in yaw coefficient 





Symbol Definition 


CLP(120) damping in roll coefficient 

CM(121) pitching moment coefficient 

CN( 122) yawing moment coefficient 

Gi 123) rolling moment coefficient 

CG1(136) center of gravity (ft) 

A(201) moment-of-inertia about missile roll axis (x-body 


axis) (slug-ft*) 

B( 202) moment~of-inertia about missile pitch axis (y-body 
axis) (slug-ft*) 

CC(203) moment-of-inertia about missile yaw axis (z-body 


axis) (slug-ft*) 


TSA(208) angle between SA and BA (rad) 

P(212) missile angular velocity about x-BA (rad/s) 

Q(216) missile angular velocity about y-BA (rad/s) 

R(220) missile angular velocity about z-BA (rad/s) 

AG(282) unit conversion lb-slug 

VXTP (286 ) missile velocity in x-TF (ft/s) 

XTP(290) missile displacement in x-TP (ft) 

VYTP( 294 ) missile velocity in y-TP (ft/s) 

YTP(298) missile displacement in y-TP (ft) 

VZTP( 302) missile velocity in z-TP (ft/s) 

ZTP( 306) missile displacement in 2-TP (ft) 

GZRO( 404 ) constant set to zero for flat earth gravitational 
field and set to ons for a spherical gravitational 
field 

ER(405) angular velocity of earth (rad/s) 

ALPO(406 ) angle between north and x-TP (rad) 





Symbol Definition 


OLAMO(407) latitude origin of tangent plane 

HO(414) distance tangent plane is from earth (ft) 
GO(415) gravitational acceleration (ft/s*) 

HREF( 501) reference altitude for midcourse guidance (ft) 
RE( 503) earth's radius (ft) 

H( 507) altitude normal to earth (ft) 

AMACH( 520) missile Mach number 

THET (521) missile pitch angle, TP (rad) 

PS1(522) missile yaw angle, TP (rad) 

PHI( 523) missile roll angle, TP (rad) 

GAMMAV( 527) vertical flight path angle, TP (rad) 

VEL( 528) magnitude of missile velocity (ft/s) 

VAT( 529) missile velocity (ft/s) 

TF( 550) program termination time (s) 

VAH( 561) computed speed of sound (ft/s) 

ALAT( 576) latitude of target position (deg) 

AZ(578) azimuth of target position (deg) 

DYNP( 581) dynamic pressure (psi) 

FLGRJ (606) constant set to zero indicates run will use 


ENGINE subroutine and when set to one run will 
use RAMJET subroutine (simplified ramjet model) 


T(932) actual time (s) 

T1(933) boost engine ignition (s) 

T2(934) commence acceleration command mode (s) 

T3(935) boost engine burn-out/port cover blow-in (s) 
T4(936) start ramjet engine (cruise) (s) 

T5(937) commence heading and altitude guidance control (s) 
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Symbol Definition 


T6(938) commence terminal dive (s) 

T7( 939) commence terminal guidance - search (s) 
T8(940) commence terminal guidance - track (s) 

DMAX (1740) maximum fin deflection (rad) 

CPP( 2669) time between printouts (s) 

ROLLO(2901 ) initial roll angle (rad) 

PITCHO( 2902) initial pitch angle (rad) 

YAWO( 2903) initial yaw angle (rad) 

STEP(2905) determines executive program flow after staging 
DOC( 2909) defines number of times COMMON will be printed 
HMIN(2911) minimum integration step size 

HMAX( 2912) maximum integration step size 

DER(1)(2913) integration step size (s) 


sla 
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I. INTRODUCTION 


Research was undertaken to convert the MOD6DF computer program re- 
ceived from NWC China Lake for use on the IBM-360 computer at tne Naval 
Postgraduate School. Once a functioning program was obtained, initial 
displacement error versus impact miss distance was investigated for the 
supersonic tactical missile. Additionally, the midcourse guidance 
reference altitude was changed to reflect a sea-skimming scenario and 
this was examined for its effect on the terminal guidance problem. 

This report not only discusses the aforementioned topics, but also de- 


scribes the missile mission requirements and the MOD6DF computer progran. 
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IZ. MISSION REQUIREMENTS 


The Supersonic Tactical Missile (STM) mission is divided into six 

phases (Ref. 1). 

* initial conditions 

* separation 

* boost 

* transition 

* cruise 

* terminal 
These divisions are based on the missile aerodynamics. The initial 
condition phase establishes the starting conditions for each launch. 
This is done while the missile is still attached to the launch platform. 
The separation phase starts when the missile is launched. The missile 
falls for approximately five seconds until the boost engine ignites. 
This initiates the boost phase which continues until the missile achieves 
Mach two. As the missile passes through Mach one, plume effects are 
encountered which the aerodynamics account for. At the end of the Ddoost 
phase, the vort covers blow in. This initiates the transition phase. 
This phase is very short and allows the debris to be ejected from inlet 
ports. The cruise phase commences when the ramjet engine ignites. This 
engine propels the missile until target impact. The terminal phase of 
the flight begins when the missile is commanded to dive from the cruise 
altitude. This phase concludes when the flight is terminated at target 
impact. 

The six STM mission phases require four phases of control. These 


control phases are: 





* separation 

* midcourse guidance 

* terminal dive 

* terminal guidance 
The four Guidance, Navigation, and Control (GN&C) system control phases 
are directly related to the mission phases. Figure i displays this 
relationship (Ref. 2). To obtain a successful flight, various types of 
control and guidance functions are required. Figure 2 (Ref. 3) illus- 
trates the guidance control mode sequencing in relation to the mission 
phases and also indicates the critical missile switching times. The 
control phases are discussed below, along with the appropriate guidance 


modes. 


A. SEPARATION 

The separation phase commences upon launch. The missile is ejected 
downward from the launch platform. Additionally, the pitch attitude of 
the missile is commanded down. Since this portion of the flight is 
unguided, the ejection force and gravity are the only forces acting on 
the missile. At launch, the missile is required to be in the attitude 
command mode. 

Five seconds into the flight, the missile pitch attitude is commanded 
up and the boost engine is ignited. The booster continues until the 
missile attains Mach two, at which time control is shifted from attitude 
to acceleration control mode. The acceleration control then requires 
the missile to maintain the normal and lateral accelerations at zero. 

The last event to occur in the separation phase is inlet port cover 
blow-in. The time delay that allows the port covers to clear is a func- 


tion of altitude. 
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B. MIDCOURSE GUIDANCE 

Midcourse guidance begins upon completion of the separation phase. 
The moment the port covers are clear, the phases shift. During this 
phase two things are required. The GN&C must guide the missile to the 
established cruise altitude and then guide the missile (in the horizontal 
plane) to the predetermined target location. 

To accomplish these requirements the GN&C system employs both 
altitude and heading control. The altitude control acts upon the pitch 
axis to drive the missile to the cruise altitude and then to maintain 
that altitude until the terminal dive phase commences. 

The actual guidance work is performed by the yaw axis. The GN&C 
system uses the heading control to steer the missile to the target posi- 
tion. To generate the steering commands a guidance law is necessary. 
The guidance law should minimize the cross-track error and the final 
lateral acceleration. A minimum cross-track error allows for mininum 
flight time and minimum displacement error when target search is initi- 
ated. Minimizing the final lateral acceleration ensures that tne seeker 
will continuously be pointing toward the target area. To insure 
accuracy, the guidance law must be able to perform these functions when 


subjected to disturbances such as wind and thrust misalignment. 


C. TERMINAL DIVE 

The terminal dive phase commences when the missile is commanded to 
dive. The actual time that this command is given is a function of the 
predetermined target coordinates. Upon diving, the missile must accel@ 
erate dowrward along a 60-degree (from horizontal) dive angle to the 
target position. Once the missile has steadied up in the proper dive 


aspect, terminal guidance commences. 
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D. TERMINAL GUIDANCE 

During the terminal guidance phase the seeker is required to search, 
acquire, and track the target. The GN&C system then uses these tracking 
signals to direct the missile to the target. Ihe seeker uses a prepro- 
gramed search pattern to locate the target. When the target has been 
acquired, the seeker then commences tracking it and also notifies the 
GN&C system that the target is being tracked. 

Target acquisition takes place when the target falls within the 
instantaneous four-degree beamwidth pattern as it traverses the scan 
pattern on the earth's surface. The computer program uses this assump- 
tion since the exact missile target acquisition mechanism has not yet 
been defined. 

Upon acquisition, the GN&C system closes the seeker tracking loop. 
The radiometer error signal is defined as the difference between the 
look vector and the line-of-sight. This signal is used to reposition 
the antenna gimbals so the look vector and the line-of-sight coincide. 
Under the closed loop operation the radiometer output is proportional 
to the line-of-sight rate and this is the signal used to guide tne 


missile to target impact. 
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III. COMPUTER SIMULATION 


The modularized six-degree-of-freedom (MOD6DF) computer program was 
developed by the Litton Systems, Inc., Guidance and Control Systems 
Division (Ref. 4) to be used in analyzing missile guidance and control. 
The program uses a building block approach, where each module corresponds 
either to a missile subsystem or an environmental system. In its 
original form, the program is used primarily for terminal guidance of 
air-to-surface missiles. The Naval Weapons Center, China Lake modified 
the original program to specifically apply to the Supersonic Tactical 
Missile (STM). The modification also allows the user the capability of 
simulating any portion of the missile flight trajectory. 

The MOD6DF computer program consists of four main decks and one 
auxiliary deck. The main decks include the executive programs, 
operational subroutines, modules, and input data deck. All the sub- 
programs utilized in the integration and the sequencing of the modules 
are contained in the executive programs. The operational subroutines 
are used by the user to control the program while it is running. The 
psysical system and the environmental subroutines are included in the 
modules. The input data deck contains all the information the executive 
programs need to execute the desired run. Lastly, the auxiliary deck 
consists of all those subprograms (subroutines and functions) that are 
required by the modules. 

Once the user starts making simulations, he must concern himself 
with program sequencing. Froper sequencing is required to ensure a valid 
run is achieved. Since it is of such importance, sequencing will be dis- 


cussed as the final part of this section. 
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A. EXHCUTIVE PROGRAMS 

The executive programs are the main core of the MOD6DF program. 
These subroutines have various functions in setting up, sequencing, 
integrating, and resetting the program. Since these functions are 
required by any type of analysis done, the user should never have to 
change any of these subroutines. Reference 5 should be consulted if the 
user desires more information about the content of any of the subroutines 
described herein. 


le Lero 





Subroutine ZERO is used to set all the COMMON(3415) locations to 
the value zero. This is done to ensure that no erroneous information is 
used in the simulation. 

Ze Cinpt i 

Subroutine OINFT1 is the basic input routine for the MOD6DF 
program. The normal input is from punched cards. However, inputs may 
also be read from tapes. The basic input cards will be discussed in the 
section covering the input data deck. 

3. Auxi 

Subroutine AUXI is used to call the initialization modules for 
each run. Additionally, it sets up the list of state variables which 
are used in AMRK. 

4, Auxsub 

Subroutine AUXSUB is used to call the dynamic modules. In 
calling the dynamic modules, AUXSUB sets and resets the lists needed for 
AMRK and the COMMON(3415) storage cells with the most recent values of 


the state variables and their derivatives. 
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5. Amrk 
Subroutine AMRK is the integration subroutine. It uses a point- 
wise first order Runge-Kutta method. All the state variables to be 
integrated must be listed in COMMON(3415) and also must appear in a 
processing list. 
6. Reset 
Subroutine RESET is used to reinitialize up to fifty input 
parameters. This is done prior to the start of any repeated runs. To 
indicate which parameters are to be reset, the number one is punched in 
columns 46-60 on the respective type three card. 
7. Return Group 
The return group is a collection of all the unused modules. 
Since all the modules and their initialization modules are called by 
AUXI and AUXSUB, the unused ones must still remain in the deck. These 
are required to ensure proper linking when the computer attempts to 
link all the subroutines. All the subroutines in this grouping contain 
three cards. The three cards are the subroutine title, and a return and 
an end card. 
8, Subli, Subl2, Subl3 
These subroutines are used to call the operational subroutines 
that are required. They call the routines in the order prescribed by 
the input data deck. The number at the end of each subroutine title 
indicates which operational subroutines it can call. For example, 


SUB11 ean call STGE1. 
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B. OPERATIONAL SUBROUTINES 
The operational subroutines provide the user with control of the 
program while it is running. The order in which the operational sub- 
routines are called is specified by the input data deck. Since these 
routines assist the user in controlling the simulation, they can be 
reprogrammed. However, it is advised that they not oe changed until the 
user has become quite familiar with the overall operation and sequencing 
of the MOD6DF program. For more in depth knowledge of the operational 
subroutines the user should consult Ref. 6, and the computer listing, 
which is contained at the end of this report. 
1. Inpti, Inpt2, Inpt3 
These subroutines are available for new inputs during the simu- 
lated flight. The only one presently used in the program is INFT1. It 
is utilized to input a namelist file which is used by ENGINE and is 
described there. 
2. Oupti, Oupt2, Ount3 
These subroutines allow for the print-out of up to fifty dif- 
ferent variables during the flight simulation. OUFT1 is not utilized 
in the deck. OUFT2 is used when the desired outvnut is to be put on 
tape. OUPT3 is the basic output routine. It prints the desired output 
on regular computer paper. 
3. Stgel, Stge2, Stze3 
These subroutines allow for proper staging, run termination, 
etcee Presently STG#i is not being used. STGE2 is being used as the 
staging initialization subroutine. STGz3 is the primary subroutine of 
this group. It stages when impact with the earth is made, when the 
final time, TF(550), is reached, or when LCONV(2672) is set equal to 


two. All the tolerances for staging are listed in STGS3, 
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ib, Chum pcmene . Cntr) 


These subroutines allow the external dynamic control inputs to 

the modules. In the MOD6DF program these are not used. 
5. Rndml, Rndm2, Rndm3 

These subroutines allow random noise to be added to the state 
variables generated in the modules. RNDM1 is not used in the program. 
RNDM2 is used as the initialization subroutine, while RNDM3 provides 
continous noise values. These subroutines produce correlated noise 
values for as many modules as required. The noise values remain fixed 
during each individual integration cycle. 


6. Auxal, Auxa2, Auxa3 
Auxbi, Auxb2, Auxb3 


Auxci, Auxc2, Auxc3 


These subroutines are auxiliary routines that allow for external 
input, output, control, etc. of the modules. At the present time, none 


are utilized in the MOD6DF program. 


C. MODULES 

The modules are of prime importance to the user since they represent 
the ‘model® of the dynamic system. In general, the model is described 
by ordinary non-linear time-varying differential equations with both 
random and deterministic forcing functions. The user must first reduce 
these equations to an equivalent system of first order equations, which 
can then be described by each module. Generally speakinz, the physical 
system is so complex that this would be impossible to do. However, due 
to the modularity of the MOD6DF program, the user can think of each 
modules as a completely independent system described by the equations 
within that medule. 


There are thirty-six possible modules divided into five functional 


24 





categories. Each group is identified by a letter which pertains to that 
groups function; A (airframe), C (computers), D (dynamics), G (geo- 
physical), S (sensors). A complete printing of each module is contained 
in the computer program listing. 
i. Airframe 
ae Subroutine Al 
This is the aerodynamic forces and moments modules. It cal- 
culates all the necessary forces and moments in body axes. These 
values are then used in the computation of the dynamics. 
be Subroutine A2 
This is the missile aerodynamic coefficient module. It cal- 
culates the required coefficients using the information stored in the 
BLOCK DATA. Using the timing inputs, this routine computes the coef- 
ficients for the different effects. ‘Some of the effects accounted for 
are; plume, separation, and control surfaces effectiveness. With this 
done the total coefficients are determined. 
ce Subroutine A3 
[This is the missile propulsion module. The timing inputs are 
used to determine whether the missile is in free fall, boost, transition, 
or cruise phase. With this determined the correct engine subroutine 
(BOOST, RAMJET, ENGINE) can be called. Three variables are calculated 
using the body axes as the frame of reference. They are the missile 
thrust is all coordinate directions, the principal moments of inertia, 
and the missile weight. 
d. Subroutine A4 
This is the fin actuator module. The four control surfaces 


commands are calculated as either ideal actuators or as second-order 
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ones. In addition to control surfaces commands, the rate of change of 
yaw, and roll are computed. 
@e Subroutine A5 
This module is part of the return group. 
2e Computers 
ae Subroutine Ci 
This is the autopilot module for the STV-G. It uses the 
cruise engine ignition time to divide the routine into boost/separation 
and cruise pnases. These two phases use different algorithms to 
calculate the turning moments in pitch/yaw and roll. 
be Subroutine C2 
This is the guidance command module. Using the timing inputs, 
this routine is divided into separation/boost, dive/climb, cruise, 
terminal dive, and terminal homing sections. Each section uses slightly 
different algorithms to calculate the guidance commands to maintain the 
proper flight profile. 
Ce Subroutines C3 =- Ci0 
These modules are part of the return zroup. 
3. Dynamics 
ae Subroutine Di 
This is the translational dynamics module. It computes the 
total acceleration in body axes and then converts them to the tangent 
plane reference. Then, accounting for aerodynamics, thrust, gravity, 
and coriolis, the velocity and acceleration are calculated. 
be Subroutine D2 
This is the rotational dynamics module. With the principal 


axes as a reference, this subroutine computes the body angular rates and 


the attitude direction cosines. 
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ec. Subroutines D3 - DS 
These modules are part of the return group. 
4, Geophysical 
ae Subroutine Gl 
This is the gravitational and coriolis acceleration module. 
It calculates the gravitational acceleration using one of two fields. 
The user specifies the field to be used by an input card. To use a 
flat-earth gravitational field, GZRO(404) must equal 0.0, and to use a 
spherical gravitational field, GZRO(404) must equal 1.0. 
b. Subroutine G2 
This module is part of the return group. 
¢c. Subroutine G3 
This is the air data module. It computes the velocity, in 
all three coordinate directions, with respect to the air mass. These 
values are then resolved into body and stability axes. This module also 
computes all the properties of air by calling subroutine AIR. These 
values are stored in their COMMON(3415) locations for use in the other 
modules. 
d. Subroutine G+ 
This module is part of the return group. 
@. Subroutine G5 
This is the coordinate conversion module. It takes the 
missile position, does a coordinat2 conversion and then it determines 
the position, velocity, and acceleration in the ECI system. 
f. Subroutine G6 


This module is part of the return group. 
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5S. Sensors 
ae Subroutine Si 
This is the homing seeker module. It simulates the missile 
seeker and computes the seeker ee and suler angle rates. Several 


flags are used to control the seeker sequencing: 


(1) FLAGS(335). FLAGS signals the start of the seeker 
search. 

(2) FLGT(317). FLGT signals when the seeker is locked-on. 

(3) FLGTS(347). #FLGTS signals the end of search. 

(4) FLGD(336). FLGD signals when the seeker has detected 
the target. 


(5) FLAGLT. FLAGLT signals when the target is outside the 
seeker field of vision. 
bo Subroutine 52 
This is the radiometer module. I1t takes the target position 
and the ATIGS target position and converts them from cody axes to the 
seeker axes. Using target position, the module tnen calculates the 
azimuth and elevation error signals. 
ec. Subroutines 53, S4& 
These modules are part of the return group. 
d. Subroutine 55 
This is the accelerometers and gyros modules. The user has 
the option of using ideal or digital accelerometers. To specify the 
type of accelerometer, the user must include the appropiate data 
statement in the subroutine. If ideal accelerometers are desired, FLGA 
must equal 0.0 and for digital accelerometers, FLGA must equal 1.0. 
e. Subroutines S6 - S10 


' These modules are part of the return group. 
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D. INPUT DATA DECK 

The input data deck provides the user with the means of specifying 
which operational subroutines and modules are to be utilized for the 
desired run. It also allows the user to set the starting conditions. 
In general, only the constants and the state variables must be given 
initial values. All quantities in COMMON not given initial values will 
be set to zero by ZERO. Im addition to the state variables, the upper 
and lower error bounds must be initialized. There are seven types of 


input cards, each indication a certain function. 


Type Function 
8) read/write tape 
1 operational subroutine to be called 
2 module to be called 
5, numerical input 
4 printed output 
5 parameter square and stm 
6 termination and random noise generator input 


A separate card is recuired for each subroutine, module, input, and 
output quantity. A sample computer printout of tne input data deck is 
contained in the computer output section. 
1. Type O Card 

Type O cards are used to indicate if the type 3 inputs are to be 
read from or written onto an auxiliary tape. These procedures can be 
used rather than reading the inputs from a deck of cards. A typical card 
is defined by punching a zero in column 2. The field that covers 
columns 5 - 20 is used by the user for any descriptive statements with 


which he wishes to identify the input. Column 21 -25 contains the right- 
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justified integer number of the tape transport to be used. The number of 
the first record to be read is punched in column 31 - 45. The last 
field is column 46 - 60 which contains the number of records to be read. 
The last two fields may use either fixed or floating-point notation. 
2. Type 1 Card 

Type 1 cards are used to specify which operational subroutines 
are called during the flignt simulation. This type card is identified 
by the number one in column 2. The second field is column 5 - 20 which 
contains any identifying information. This information is printed out 
when the data deck is read and allows the user to read exactly which 
subroutines were called. Column 21 - 25 contains the right-justifiead 
integer number which is the subroutine identifying number. [he orera- 


tional subroutine numbers are; 


Subroutine Subroutine Number 
UNE T ics ENPI? , iNPIT3 2 
OUPTI, OUPT2, OUFT3 3 
STGE1, STGE2, STGE3 4 


CNTR1i, CNTR2, CNTR3 
RNDM1, RNDM2, RNDM3 
AUXAL, AUXA2, AUXA3 


AUXB1, AUXB2, AUXB3 


Oo Oo NY BW WN 


AUXC1, AUXC2, AUXC3 
It should be noted that all or any of the subroutines listed under one 
one number can be called by including only one card. The cards are 
placed in the data deck in the order in which they will be called. 
This order or sequencing will be explained further in section F, 


Sequencing. A typical type 2 card is identified by the number two in 
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column 2. In general, column 5 ~ 20 should contain the module title, 


but any pertinent information is allowed. The module number is punched 
in column 21 ~- 25 and it must be right-justified. A listing of the 


module numbers follows: 


Module Module Number Module Module Number 
A1,Ail 2 D4, D4I 20 
A2,A2I 3 D5,D51 Zi 
A3,A31 4 Gi,Gil 22 
AGL, AUT 5 G2,G2I 23 
A5,A51 6 G3,G31 24 
C1,Cil C G4 ,G4I 25 
C2,C2I 8 G5,G51 26 
C3,C31 9 G6,G6I1 27 
CH, CUI 10 $1,SiL 28 
C5,65. aba S2,0cL 29 
C6,C61 12 S3o3L 30 
C7675 3 S4 ,S4I 53 
C8,C8I 14 55,051 32 
ES5CSI 15 S6,S6I 33 
Ci0,Ciol 16 S7,97L ce 
Di, DiI 17 S8,S81 35 
D2,D21 18 59,991 36 
D3,D31 19 910,5101 37 


Note that either the module, the initialization module, or both may be 
called by including only one card in the deck. A sample of a typical 


type 2 card is shown in Figure 3 (Ref. 7). 
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TYPICAL TYPE 3 CARD 


FIGURE 4 


32 





4, Type 3 Card 

Type 3 cards are used to set any COMMON(3415) location to any 
value other than zero. In general, four items must be initialized. The 
state variable initial values and any constants used in the flight simu- 
lation are the most obvious. Additionally, there are some constants 
associated with the executive programs and operational subroutines and 
the state variable upper and lower bounds which must be initialized. As 
with all cards, column 2 defines the type card and it must contain a 
three. Column 5 - 20 holds the statement describing the input. The 
user should be specific here since it will save him having to remember 
every COMMON(3415) location. The only other means of input identifica- 
tion is by column 21 - 25. These columns contain the right-justified 
COMMON location of the input. Columns 31 - 45 hold the actual numerical 
value of the input. The last field, column 46 - 60, contains the reset 
flag. If the reset flag equals one, the COMMON(3415) location and the 
numerical value are placed in the reset list. This list may contain up 
to fifty different variables. This, in the case of multiple runs, 
allows the variables to be reset to its initial value prior to each run 
without additional input cards. The sample card in Figure 4 (Ref. 3) 
shows that either fixed or floating-point notation may be used to input 
the numerical value and the reset flag. These cards need not be inputed 
in any specific order, but for ease of checkout, it is advised to place 
them sequencially by their COMMON location. 

5. Type 4 Card 

The MOD6DF program can printout a maximum of fifty variables for 

each simulation. Type 4 cards are used to specify which variables are 


to be printed. Column 2 must contain the number four to indicate a type 
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4) card. When the results are printed out headings are included. These 
headings are designated in column 9 - 20. The exact alphanumeric title 
punched will be printed at the top of each page, this need not be the 
fortran symbol used within the program. The COMMCN(3415) location of the 
output variable is contained in column 21 - 25 and must be right- 
justified. 
6. Type 5 Card 

Type 5 cards are used to indicate which variables are to be root- 
mean=-squared. These cards are similar to type i and type 2 cards. 
Column 2 must contain the number five. Any pertinent information about 
the variable to be operated on is punched in column 5 - 20. The COMMON 
(3415) location of the variable must be right-justified in columns 21 - 
25. The last field, column 31 - 45, indicates whether the root-mean- 
Square operation will occur along the trajectory or at the end. 

7. Type 6 Card 

The type 6 card has two purposes. Its first function is trivial, 
but required. The number six is typed in column 2 and the rest is 
blank. This indicates to the computer program that there are no more 
input cards. Its second function involves random inputs. This card is 
used to indicate the number of random process (noise) generator cards 
that are to be read before the input process is terminated. Column 2 
has the same information as before. Any pertinent information is con- 
tained in column 5 - 20. Columns 21 =- 25 must be right-justified and 


they contain the number of random generator cards to be read. 


E. AUXILIARY DECK 
The auxiliary deck is a collection of subroutines and functions 


required by the modules. These routines are general in nature since 
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they can be called by several modules. They are used to calculate such 
things as the properties of air, engine performance, various ratios, and 
to locate values in the many data tables. A brief discussion of the most 
common subroutines is presented. 
1. Boost 
This subroutine calculates the thrust coefficient (CT) and the 
fuel flow rate (FF) during the boost phase. These values are returned 
to module A3 and used to calculate the missile thrust in the body axes 
and the missile weight. 
2. Ramjet 
This subroutine is used during the midcourse, cruise, phase of 
the missile flight. It uses a simplified ramjet model to calculate the 
thrust coefficient and fuel flow rate. This routine is not automatically 
called. The user must designate that he wishes to use it by inputing a 
type 3 card setting FLGRJ(606) equal to one. This then sets up the 
proper stepping in module A3. 
3. Engine 
This is the primary routine during the cruise phase. It is one 
of many routines within the NWC air-breathing propulsion package. This 
entire package is utilized to calculate the engine parameters from inlet 
to exhaust. Again, thrust coefficient and fuel flow rate are eventually 
computed and returned to module A3. The user does not need to supply any 
spacial input cards to use this subroutine. If no initial value is 
inputed for FLGRJ(606) it is automatically set to zero, which indicates 
this routine is to be used. 
4, Air 


All the properties associated with the air are calculated in 
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this subroutine. This includes the computation of the speed of sound and 
the dynamic pressure. Since the missile does vary in altitude, the 
routine takes this into account as well as the latituds. Once these 
values are calculated, they are returned to module G3, 
5. Block Data 
This routine contains data tables. These tables cover para- 
meters from aerodynamic coefficients to thermal properties. The total 
data package covers specific bands within tne missile operating envelope. 
Data is stored in matrices which includes one, two, and three dimensional 
ones. These data tables are readily available and there are routines 
designed to retrieve this information quickly. 
6. Serch 
This subroutine, along with several functions, is used to 
retrieve information from BLOCK DATA. If the present overating point of 
the missile is not within one of the bands of information, then the 
tables are interpolated. The functions THRZDL, STDLI, STDLIA, and TA3 do 
the interpolation of the tables. Since the tables are of various 


lengths, these functions are very general. 


F, SEQUENCING 

Sequencing is very important in the running of the MOD6DF program. 
Care must be taken to ensure that the modules are processed in the cor- 
rect order at each step. This is essential to eliminate the use of 
obsolete values from the last cycle. An exception to this vroblem is 
the state variables. These are updated simultaneously by the integra~ 
tion algorithm. Any module is capable of using the most current value 
of these, no matter what the order of processing. 


To help remedy this problem, module diagrams were devised. Module 
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diagrams aid the user in maintaining the proper flow of variables into 
and out of the modules. To design a module diagram start with a box. 
This box will contain all the equations for a specific module. The 
standard procedure for showing inputs and outputs is to use arrows 
pointing in or out from either the left or the right side. The example 
in Figure 5 (Ref. 9) shows this technique. The arrows pointing in from 
the left indicate variable inputs from other modules. Arrows pointing 
out to the right indicate output going to either other modules or 
program output. The last set of arrows points in from the right. These 
show the constants brought in directly or indirectly from the initializa- 
tion module. Since each arrow represents a variable, thsy must be 
defined. The usual means of labeling the arrows is to use the variable 
fortran symbol and in parenthesis its COMMON(3415) location. 

Hach variable usually has only one COMMON location associated with 
it. In the case of ‘state’ variables this is not true. State variables 
are defined by four consecutive COMMON locations. The first is for the 
derivative of the variable. The second and the third locations hold the 
ever and upper bounds, respectively, of the integration error. The 
last one contains the variable itself. 

Once all the required module diagrams have been completed, they can 
be combined to get the overall processing order. The usual processing 
order, shown in Figure 6 (Ref. 10) is to start with the geophysical 
group, then proceed to sensors, computers, airframe, and finally 
dynamics. Within each group is a usual processing order and this too is 
shown in Figure 6. This process for determining the program sequence 


will eliminate the use of any obsolete values in the computer run. 
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IV. COMPUTER PROGRAM CHaCKOUT 


The basis for the research was the MOD6DF computer program from NWC 
China Lake. The total package received from them consisted of a listing 
of the program and an uninterpreted deck of cards. The initial step was 
to input the cards in small groups into the computer and then to examine 
the source listing. This listing revealed that the original deck was 
punched in BCD. This fact was easy to determine since several characters 
were changed (Ref. 11). The library routine NeiwDfK was used to translate 
the cards from BCD to EBCDID. 

Once the translation was completed, an attempt was made to compile 
the new deck. This produced an output which contained many syntax errors, 
Misse errors were divided into two major groups. Une effected the use of 
quotation marks in FORMAT and comment statements. The other one, the 
more difficult, effected the DATA statements in the BLOCK DATA subroutine. 

The problem with the DATA statements was due primarily to tne differ- 
ence in the compilers used. The compiler at NWC was much newer and 
allowed for the use of more sophisticated inputs. The compiler at NFS 
only allows a data set to start with the first element. This required 
the rewriting of many data groups. To complicate these revisions, a 
limit of nineteen continuation cards is also imposed. These restrictions 
demanded not only the rewriting of many data sets, but also the forma- 
tion of two new ones. 

With the corrections finally completed, the computer would then 
compile the program. The next step was to link all the subroutines 
together. The first attempt was unsuccessful. Inadvertently, the 
subroutine INTR20 had been omitted from the original deck. Using the 


program listing, the contents of INTR20 were typed and included in the 
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main deck. With this addition the program would now compile and link. 
Now that the program would compile and link, attention was turned to 

getting a good simulation. To facilitate this process, two sample 

outputs were obtained from NWC. These outputs cover one missile flight 

which is broken down into a midcourse guidance and a terminal guidance 

s imulation. Hereafter these will be referred to as midcourse baseline 

and terminal baseline, respectively. Using the initial conditions from 


the baseline models, it was hoped that the outputs could be duplicated. 


A. MIDCOURSE GUIDANCE FLIGHT 

The initial run, using the midcourse baseline inputs, revealed an 
overflow problem with the dynamic pressure ( QD(508) ). Us ing the 
traceback procedures outlined in the Users Manual from the W. R. Church 
Computer Center, the problem was confined to subroutine AIR. The 
problem turned out to be a translation error. The symbol PF? (J - zero) 
had translated to PG? in one place and PO in another. This error caused 
the program to use a value left in that memory location from a previous 
run to calculate the dynamic pressure. With this problem remedied, tne 
program could progress a little farther. The next stumdling block 
appeared as a divide check. These errors were resolved oy introducing 
patches that would bypass a statement that tried to divide by zero. 
Once bypassed, that quantity would be set equal to zero. This was the 
normal procedure of the computer, but it would stop the run after ten 
such errors. Having corrected all these errors, output was obtained 
which covered the desired seventy seconds of flight. 

When the output was examined a major switching problem within sub- 
routine A3 (missile propulsion module) was found. The midcourse 


baseline utilized a simplified ramjet model, but the output was not. 


44 





(99S )ANLL 


09 


OS 


££ daunole 


(ASMNOOCIN) ANIL *SA AGNLIVIV 


07 Ob 


OT 


OT 


Od 


O€ 


07 


(13%) AANLILIV 


42 








Rather than use RAMJET, the output disclosed that KIGINE had been used. 
The problem was in the logic statement, FLGRJ 0.0. The program was not 
making the desired switch to the simplified ramjet model. To remedy the 
problem, the data card FLGRJ(606) 1.0, was added to the input data deck. 
Another problem was discovered which occured between 14.0 and 14.5 
seconds. During that time period the missile angle of attack ( ALPHA 
(330) ) exceeded ten degreese When ALFHA exceeds this angle the engine 
is turned off. With the engine off, no thrust is produced causing the 
forward velocity ( VATP(286) ) to decrease. This limitation was removed 
from the program. After eliminating these proolems a flight trajectory, 
Figure 7, was obtained which closely resembled the output of the mid- 
course baseline. A random sampling of the outputed variables were 
compared for exactness. The differences noted were due primarily to 
computer round-off error. 

To further verify the accuracy of the two simulations, three param- 
eter ( TXBA(073), VXTP(286),and YTP(298) )were chosen as representative 
values for the runs. To obtain a feel for the run, these parameters 
are plotted in Figures 8, 9, 10. Additionally, random time samples are 
tabulated in Tables I, lI, Iil. 

The thrust ( TXBA(073) ) plot, Figure 8, can easily be divided into 
four mission phases. During the separation phase (0.0 - 4.5 seconds) 
the thrust is zero. This is expected since neither of the engines have 
ignited. This is followed by a rapid increase in the thrust. The 
maximum thrust, 27800 lbs, happens during the boost phase (4.5 = 9.5 
seconds). At 9.5 seconds the boost motor stops and the transition phase 
occurs for the next 0.3 seconds. During this time the thrust decreases 


rapidly since neither engine is on. Once the engine port covers are 
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clear, the ramjet sustainer ignites. This action initiates the midcourse 
phase (9.8 - 70 seconds) and the thrust remains relatively constant at 
1200 lbs. 

The missile forward velocity ( VXTP(286) ) can be directly related to 
the thrust. During separation no thrust is produced, tnerefore the 
velocity is zero. As the missile is boosted to Mach two, the velocity 
increases and approaches its maximum value (2800 ft/s) just prior to 
boost engine shut-down. As the gustainer engine port covers clear, the 
velocity decreases slightly. Once the ramjet engine ignites the velocity 
profile is displayed in Figure 9 and it closely resembles that for the 
midcourse baseline. 

The last parameter, y-displacement ( YTP(298) ), was included to 
demonstrate the accuracy of the guidance system Figure 10 shows that 
prior to attaining the cruise altitude, the missile wanders off track. 
Once the guidance system. is activated, it makes the necessary corrections 
to return the missile to the planned flight path. Table III shows that 
the y-displacement corresponds to tne baseline case and at the conclusion 
of the run the off-track error is down to 12.03 feet. 

All the information obtained from the new run was checked against 
the midcourse baseline. The results showed that the two simulations are 
within acceptable limits. With this milestone completed, the program 


checkout could proceed to the terminal guidance flight. 


B. TERMINAL GUIDANCE FLIGHT 

| To checkout the terminal guidance portion of the MOD6DF program, the 
initial conditions from the terminal baseline were used. This required 
changing about a dozen input cards in the midcourse input data deck. 


The first run attempted turned out successful. This was due to the 
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elimination of all the errors during the midcourse checkout. The result- 
ing flight trajectory, Figure 11, is a very good illustration of the 
desired profile. To check the accuracy of the simulation, the same 

three parameters were chosen. Tables IV and V show that this run is 
almost identical to the terminal baseline. The difference is primarily 
due to computer round<-off error. Figures 12 and 13 show the overall 
flight. performance of the thrust and forward velocity, respectively. 

For this simulation, the y-displacement ( YTP(298) ) does not 
compare favorably with the terminal baseline. Figure 14 shows that be- 
tween ten and twenty seconds in the flight, ITP develops rapid 
oscillations. By examining Table VI, it is noted that during this 
interval the two simulations differ the most. No reason for this 
discrepancy could be found. These oscillations do effect the flight by 
increasing the time of flight from 25.825 to 29.09 seconds and by 
increasing the miss distanre. The imposed time constraints force this 
problem to be overlooked and to direct attention to the accuracy of the 
impact. Since the impact is within 0.5 feet of the target location, it 
is concluded that the terminal guidance simulation works properly. 

With both segments of the MOD6DF computer program working correctly, 
thoughts could now turn towards the actual experimental runs. The 


results of these runs are discussed in the next section, 
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V. TERMINAL FLIGHT DISTURBANCES 


With the checkout of the MOD6DF computer program completed, many 
ideas were discussed concerning alterations to the simulations. The 
three modifications decided upon were: 

* examine target miss distance when changes were made to the 
initial x-y-z conditions, 

* examine the terminal flight profile, range, and miss distance 
when the cruise altitude was reduced to approximate a sea skimming mode, 

* examine target miss distance when random noise is applied to 
the missile homing seeker. 

To understand the changes and results, one must be familiar with the 
frame of reference. The origin of the reference frame travels from the 
launch platform to the target location at the cruise altitude. This is 
called the tangent plane reference system. Displacements in the 
x-direction ( XTP(290) ) are measured from the launch platform in the 
direction of the target. Y-displacements ( YTP(298) ) are measured left 
or right of the ideal flight path, in the tangent plane. Any vertical 
displacement ( ZTP (306) ) is measured normal to the ideal ite. path, 
with down being the positive direction. Using this frame of reference 
the information in Table VII is easier to understand. 

The first modifications demonstrate the effect of chanzing YTP. 

YTP was increased until a result was reached that was unsatisfactory. 
From the results in Table VIII, the maximum value of YIP was determined 
to be 1800 feet. This conclusion corresponds with the results in 
reference 13. Since the missile is symmetric, it was also concluded 
that moving YTP either right or left would give the same results. 


The fifth run simulated a drop in the missile's altitude. ZTF was 
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jinputed as 2000.0 feet. As the eee progressed, the missile 
developed the necessary commands to climb and regain the desired cruise 
altitude (Figure 15). It then commenced its terminal dive. The target 
was detected, acquired, and impact followed. This initial condition 
modification produced the largest miss distance (0.457 ft.). 

Figure 16 was presented to show the relative locations of the miss 
distances for each change. Table VIII reveals that up until 1300 feet, 
changes to YIP created very little variation in the miss distance. It 
should be noted that the larger the displacement became, the greater the 
time until impact. The increase in time was necessary to allow the 
missile to acquire the target and then compute the required acuator 
commands to impact the target. This concluded the investigation of the 
effects of initial displacement error on target miss distance. 

The idea behind changing the cruise altitude was to try and simulate 
a sea skimming profile. RAMJET contains tables that allow for four 
cruise altitudes. Of these four altitudes, the lowest (500 ft) was 
chosen even though it is high for a sea-skimmer. To run this simulation 
only two input data card changes were required. HO(414) and EREF(501) 
had to be set equal to the desired altitude. The resultant flignt path 
is shown in Figure 17. 

The simulation produced an error-free output, but at first glance 
the results appeared unacceptable. Still trying to compare miss dis- 
tances, the downrange distance (XTP) was fourd to be 197916 feet. 
Comparing this to the terminal baseline (234999) resulted in an extreme 
error. It was then realized that the missile must expend more fuel at 
this altitude to attain the same speed. Therefore, the results could be 
correct. As further verification of the correctness of the run, it was 


assumed that the missile weight at impact should be fairly equal. The 
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simulation impact weight was 1126.34 lbs while the terminal baseline 
weighed 1113.12 lbs. Since no conflicting information could be found 
in reference i2, it was concluded that this simulation was correct. 

The final flight disturbance investigated was the addition of random 
noise to the missile homing seeker. This simulation proved -to be unsuc= 
cessful due to invalid input format. The information cescribing the 
process is brief and difficult to understand. Further information from 
NWC China Lake is required:«to be able to successfully run simulations 
with random noise. 

With the completion of these modifications, many unanswered problems 
and questions still exist. The two most critical problems concern the 
input of random noise generators and the missile exceeding the angle of 
attack limitation when using ENGINe for cruise propulsion. However, the 
simulations did show that the MCD6DF program runs correctly at alternate 


cruise altitudes and with initial displacement errors. 
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VI. CONCLUSION 


The MOD6DF computer program from NWC China Lake was converted to 
operate on the IBM-360 computer at the Naval rostgraduate School. The 
program would only function properly when using the simplified ramjet 
model. When this model was not used, the missile angle of attack 
exceeded the maximum limit of ten degrees. This error caused the ramjet 
engine to flame out. 

Target impact errors for the terminal guidance problem were 
investigated when the initial displacement was modified. These modiifi- 
cations demonstrate the missile’s accuracy when removed from the ideal 
flight path. The results also point out that if the target falls within 
the seeker search pattern, target impact will inevitably happen. 

Random noise generation is possible with the MOD6DF program. 
However, more information discussing the parameters required to develop 
the noise is necessary. Additionally, sample noise inputs should be 
obtained that reflect the alteration to the noise subroutines. 

This research involved the preliminary investigation of the MUD6DF 
program. The program has many possible areas, concerning guidance and 
control of tactical missiles, which could be developed for future study. 
These areas not only include the unresolved problems encountered during 
this research. Additionally, reference 14 and 15 contain many examples 


of possible advanced guidance concepts. 
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